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Abstract 

It has recently been shown that in single field slow-roll inflation the total volume cannot 
grow by a factor larger than without becoming infinite. The bound is saturated exactly 

at the phase transition to eternal inflation where the probability to produce infinite volume 
becomes non zero. We show that the bound holds sharply also in any space-time dimensions, 
when arbitrary higher-dimensional operators are included and in the multi-field inflationary 
case. The relation with the entropy of de Sitter and the universality of the bound strengthen 
the case for a deeper holographic interpretation. As a spin-off we provide the formalism 
to compute the probability distribution of the volume after inflation for generic multi-field 
models, which might help to address questions about the population of vacua of the landscape 
during slow-roll inflation. 



1 Introduction and Conclusions 



The future evolution of our Universe appears to be dominated by a phase of accelerated expan- 
sion [1] |2]. The data from the cosmic microwave background strongly suggest that also in the 
distant past our Universe experienced a phase of slow- roll inflation |3l HI |5]. Furthermore, in- 
dependently from the actual observations, the study of the physics in a accelerating universe is 
interesting in its own, as still several theoretical aspects are poorly understood. 

A particularly interesting set-up is represented by slow-roll eternal inflation [6l [TJ |8]. In this 
case the scalar potential is so fiat that quantum fluctuations dominate over the classical rolling 
of the scalar field. In this limit the scalar field becomes free, so despite quantum fluctuations 
dominate over the classical evolution the system becomes exactly solvable [9]. In this case there 
is a finite probability to generate infiaton trajectories going uphill the scalar potential, which can 
make inflation last forever. Quantum fluctuations are able to completely change the future causal 
structure of space-time, similarly to what happens in false vacuum eternal inflation [IOl[Il] or in the 
presence of a black-hole [12] . Indeed, these are the only three known solutions of general relativity 
where quantum effects induce such a spectacular behaviour. In particular, in the eternal inflation 
case, the result is even more dramatic because the background geometry becomes completely 
stochastic. 

The dependence of the future causal structure of space-time on the dynamics of the theory is 
a delicate issue already at the semiclassical level. The problem is even more acute in full quantum 
gravity, where local observables are not well defined and its proper formulation in known cases 
(such as the string theory S-matrix and AdS /CFT) strongly relies on well defined asymptotic 
boundaries. To make the problem even more involved, unitarity, and in this particular context, 
holography also seem to suggest that space-time regions causally disconnected by an horizon are 
redundant, complementary. 

These problems are particularly relevant in the framework of the landscape. If, as suggested by 
string theory, quantum gravity possesses a landscape of (meta-stable) vacua, our Universe may be 
doomed to deal with eternal inflation. While vacuum tunnelling from meta-stable de Sitter vacua 
seems the natural mechanism for generating eternal inflation in the landscape, slow-roll eternal 
inflation represents a unique framework to reliably study the onset of the transition between the 
non-eternal and the eternal inflating regime. In these type of models indeed, it exists a controllable 
parameter (the flatness of the potential) that smoothly interpolates between the two phases. It 
basically allows the study of the de Sitter phase with a tunable parameter. 

In this context a bound was found in [13] for any model of inflation in the non-eternal phase: 
the number of e-folding is always bounded by the de Sitter entropy Sds at the end of inflation. 
The relation with the entropy seems to suggest a connection with the holographic bounds. It is 
indeed in close analogy with the bound on the validity of the effective field theory description 
of the Hawking radiation quanta in the black-hole evaporation process: after a time of order the 
black-hole entropy in units of the curvature scale the EFT predictions for the Hawking entropy 
start deviating by order one from what expected from the unitary evolution — information starts 
coming out and complementarity emerges. A similar time-scale shows up in slow-roll inflation 
exactly when the eternal inflation phase is approached. 
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It was found later in [9] that in single field slow-roll inflation the phase transition to eternal 
inflation is sharp: at the critical value of the parameter 

3 

a finite probability of creating an infinite volume after inflation develops. This suggested the 
existence of a sharp bound also for the number of e-foldings. Because of quantum fluctuation 
however, the number of e-foldings is not a well defined object, as it can fluctuate from point 
to point in space. The invariant quantity is the total volume of the Universe at the reheating 
surface, which matches e^^ in the absence of quantum fluctuations. Quantum fluctuation makes 
this quantity stochastic too. A probability density for the volume V can nonetheless be defined 
and actually computed. This was accomplished in [13] where a sharp formulation for the bound 
was also found, namely: 

The probability of producing a finite volume V larger than e^isl"^ vanishes up to non-perturbative 
quantum gravity effects. 

Notice that the bound applies not only in the non-eternal regime {Vl > 1), where the volume 
is always finite, but also in the eternal regime, where there is still a non-zero probability to have 
a finite volume. 

The existence of a sharp bound and the connection with the entropy clearly cries for an 
holographic interpretation. Is this bound really set by complementarity? In that case is there a 
meaning of the factor | appearing in the bound? We are not able to give a definitive answer to 
these questions yet. It seems however that if the answers to these questions are positive then the 
bound should be universal, including the factor |, as, for example, in the Bekenstein-Hawking 
formula. 

In this paper we test the universality of the bound against three different generalizations of 
the models studied in [H] by changing the number of space-time dimensions (from 4 to D), by 
taking into account the effects of higher-dimensional operators in the action, and by considering 
more inflaton fields. 

In the first test the dependence on the number of dimensions appears in a non-trivial way in 
a number of quantities entering the calculation of the bound, such as the relation between the 
entropy and the horizon area, the Friedmann equations, the quantum fluctuations of a inflaton 
field in de Sitter space, etc. We will show that, surprisingly, when the bound is written in terms 
of the total volume and the cISd entropy the dependence on the number of dimensions cancel out 
leaving the bound unchanged, with the same factor as in four dimensions. 

The second test we performed is with respect to higher derivative terms in the action. These 
terms not only change the equations of motion for the metric and the inflaton, thus the Friedmann 
equations, but also the expression for the entropy, which is not given anymore by just the horizon 
area in Planck units, and the size of the quantum fluctuations of the inflaton. Remarkably we find 
that, independently of the corrections considered in the action, the bound is not affected, included 
the factor |, which stays universal. 

The third test requires more effort. We have to extend the formalism developed in [H] to the 
case of multifield inflation. In particular the possibility to have extended regions in field space 
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where inflation can end triggers difficulties both at the technical and the conceptual level. We will 
show that the main formula that in the single-field case was a simple Laplace transform of the 
solution of a non-linear differential equation maps in this case into a functional integral transform 
of the solution of a non-linear partial differential equation. These more involved formulae allow 
us to calculate not only the probability distribution of the total volume after inflation, but the 
multi-variate probability distribution of each different kind of volume associated with the different 
reheating points in field space. The possibility of exiting inflation in different places in field space, 
where in particular the Hubble scale and thus the entropy are different, forces us to generalize 
the definition of the bound too. A conservative assumption is to assume that the total volume 
be bounded by the largest possible entropy on the reheating surface. We believe however that a 
stronger version of the bound actually holds. 

In all the realizations with finite total volume, the probability of producing a particular volume 
with values of the inflaton field within a given region X of the reheating surface and larger than 
the largest value of e^"^^^"^ on the same region X vanishes up to non-perturbative quantum gravity 
effects. 

In particular given any reheating region in field space the volume produced after inflation with 
that kind of inflaton values will be bounded by the corresponding de Sitter entropy. 

Because of the complexity of the formulae involved we are not able to test the validity of the 
bound in the most general case, however we will discuss the two simplest examples where we are 
able to take the calculation till the end. The two examples correspond to two-field inflationary 
models with a constant slope of the potential and where the reheating region is a straight line 
normal or at angle with respect to the slope direction. We will be able to calculate the average 
volume distribution as a function of the reheating point and show that the bound is indeed 
satisfied. Interestingly the factor ^ is still universal and the bound is saturated only when the 
slope is actually orthogonal to the reheating surface. This configuration correspond to single field 
inflation with a spectator field. This result suggests that the presence of extra fields only makes 
the bounds stronger, by drifting the inflaton trajectories towards region of higher entropy where 
the bound is less constraining. 

It would be nice to test the bound also for more complex field configuration, but we are 
limited by our ability to solve the corresponding equations. Another interesting effect that we 
have neglected so far is the inclusion of slow-roll corrections. While they are expected to be 
very small, after all in the eternal inflation limit the slow-roll parameter is tiny, when considering 
inflaton trajectories that start saturating the bound, the integrated slow-roll effect may be non- 
negligible, the tiny slow-roll parameter being compensated by the large field excursion. Arguments 
were given in [T3] why this effect should not alter the nature of the bound, still a dedicated study 
would be worthwhile. 

Finally we would like to mention a suggestive relation between our bound and some thermo- 
dynamic relations recently found for systems out of equilibrium (see e.g. [15]). The variation of 
the entropy in a closed thermodynamical transformation can be thought of as a measure of the 
irreversibility of the process. If we consider the cyclic transformation of a piston moving up and 
down, the process will be the more irreversible the quicker we move the piston. Analogously in 
eternal inflation, we can think of de Sitter space as an equilibrium state, and inflation as the irre- 
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versible process associated to moving the scalar field along the potential. The slower the infiaton 
moves, the more inflation comes closer to becoming reversible and also closer to de Sitter space. 
Let us indeed consider the equality AS^s = 12QNc with Nc the classical number of e-foldings, 
that is described later in the text. This equality is usually taken to represent the second law of 
thermodynamics AS = i^Q/T)^^^, with (SQ/T)^^^ = 12flNc [HI [H]. This interpretation is a 
bit puzzling to us, because it may look like that all classical solutions of General Relativity are 
reversible. On the other hand, the presence of an horizon suggests the presence of some irre- 
versibility. In particular, slow roll inflation is irreversible at the classical level, and the flatter we 
make the potential, the closer we should get to reversibility. An interpretation of our formulas 
that would be more coherent with this intuition would be to use the bound > 1 for the non- 
eternally- inflating potentials, to say that whenever we do not have eternal inflation, AS > 12Nc, 
with the inequality being saturated at the phase transition to eternal inflation = 1. It would 
be very interesting to establish such a connection between slow-roll inflation and thermodynamics 
on solid grounds, for it may give deeper understanding about de Sitter space and the onset of the 
eternal regime. 



2 Mini-review of known results 

It has been shown in [15] that many information about the phase transition from non-eternal 
to eternal slow-roll inflation are encoded in a rather simple formula. We briefly review here the 
results of that paper, which we refer to for details. 

In fact the probability distribution p{V, r) for the volume of the reheating surface V, with 
(p = Ht/{2\/Q-k) the starting value for the inflaton fleld, is given by 



P{V,T) = —1 dzfir,z)e^'' , (1) 

2m 7o+-ioo 



where at leading order in the slow-roll approximation, / is the solution of the following differential 
equation (see also [IB]) 

d'Jir; z) - 2Vndrf{T; z) + /(r; z) log[/(r; ;.)] = , (2) 

with 

^--h^ (3) 

and boundary conditions 

f{Q-z)=e'% (4) 

drf{r-z)\r=r,=Q. (5) 

The flrst boundary condition corresponds to the end of inflation at r = = 0, while the second 
condition at r = is a barrier condition, to make flnite the allowed fleld space for the inflatoiiEI. 



^Thc barrier point T}, can be thought of as the region where the energy density of the inflaton potential becomes 
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The meaning of the differential equation is more manifest when rewritten in terms of (j) as 
follows: 

This is a modified Fokker-Planck equation: the first term is the normal dispersion term due to 
the quantum fluctuations of the inflaton field in de Sitter space (A0^/At = H^/ (47r^)); the second 
term is the drift induced by the tilt of the scalar potential; the last term encodes the volume 
growth from the de Sitter expansion. 

When rewritten in terms of r, the differential equation ([2]) only depends on the single dimen- 
sionless parameter Q — a combination of the rate of quantum fluctuations (A0^/At), the classical 
rolling (0) and the Hubble expansion {3H) — which controls the different phases of slow-roll infla- 
tion. 

Despite an analytic expression for the solution to eqs. ([T]) and for piV; r) is not available, 
the behaviors for the different regions of the parameters V and Q can be derived. Moreover all 
the moments, given by 

(0 = (-ir5:/(r;^)|.=o, (7) 

can be computed analytically, as eq. ([2]) becomes hnear for these quantities. 

For example the expression for the average volume in terms of the classical number of e- foldings 

Nc = H(p/\^\ reads 

{V) = e^^-^> = e'^^^i+Vi-vo ^ > 1. (8) 



In the classical limit, Q — t- 00, quantum fluctuation become irrelevant, and p{V] r) approaches 
a delta-function picked around the classical value Vc = e^'^". When f2 — )■ 1 quantum fluctuations 
become of order one and the average volume gets large corrections, increasing to {V) = e^'^". At 
the critical value f2 = 1, the phase transition to eternal inflation occurs and the average volume 
starts diverging. 



3 Universality of the Bound in i^-dimensions 

The first test of the eternal inflation bound we present here is its universality on the number 
of space-time dimensions. We remind that the bound follows from two ingredients: a classical 
one, which determines the largest non-eternally inflating classical trajectory, and a quantum one, 
which compute the deformation due to quantum fluctuations. As described in the previous section 
the second computation boils down to solve a differential equation and performing a Laplace 
transformation. The transition to eternal inflation happens when the only parameter of the 
differential equation Q crosses 1. The expression for the differential equation given in eq. (E]) 
is also valid in D space-time dimensions except for the coefficient of the last term, the Hubble 
expansion coefficient, which now becomes {D — 1)H. The expression for A0^/At in D-dimensions 

Planckian, or where the potential becoraes very steep, or the fixed point of a symmetric potential, such as the peak 
in the top-of-the-hill inflationary model. The limit of arbitrary far barrier, Tb 00, can be subtle and special care 
must be put in doing such a limit, see [21 [Tl] for details. 
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is also different, it can be extracted from the coefficient of tlie linear term in the 2-point function 
of a scalar field in dS^, 

TjD-l 

7rfc)(z)_i) 



l2 



where 

- ^ 

is the d-dimensional solid angle. The differential equation can be easily brought back to the 
form 02]), by defining 

o _ At _ 719(^-1) 0^ 

2{D - l)H A02 2{D - 1) i/^ ' ^ ^ 

r =cl>y/2{D-l)H^j^ = 2{D - l)VnN, . (11) 

Since the differential equation governing quantum fluctuations has the same form in D di- 
mensions, the same will be true for the solution when expressed in term of r. In particular the 
expression for the average volume will be 

{D-l)N^ 



{V) = e i+v^T^ , (12) 

with f2 as defined in eq. (fTOj) . Again for i7 — )■ oo the classical result is recovered {V) = e^^^^^^", 
while the value at the eternal inflation transition (at f2 = 1 as eq. (E]) is formally unchanged), is 
(V) = ^^iD-i)Nc ^ again quantum fluctuations increase the effective number of e-foldings up to a 
factor 2 with respect to the classical value when the phase transition is approached. 

In order to test the bound we need now to work out the relation between the classical number 
of e-folding and the de Sitter entropy in Z)-dimensions. The difference of de Sitter entropy between 
the start and the end of inflation can be written as 



Send rJ^end 1 

Sstart ^ ^start 



AS./ dS^I l^dA^-r^R^^^dK. (13) 



where we used 

for the de Sitter horizon area. From Friedmann equations in D dimensions we also have 

H = -^^\ (14) 
D-2^ ' ^ ' 

which allows us to write 



AS= I 27fQ^D-i) j^dN'^. (15) 
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We can now use eq. ffTOl) to rewrite in terms of Q and we finally get 

AS = A{D - l)ndNc = A{D - l)nN'^ , (16) 
Jo 

S,nd>^{D-l)QN,, (17) 

where we used the fact that at leading order in the slow- roll parameter Q =const. The result 
above gives a bound on the number of classical e-foldings allowed in non-eternal inflation when 
f2 > 1. The maximum number is achieved at the phase transition Q = 1, where Nc < Send/ (4(-D — 
1)). Quite non-trivially the complicated dependence on the number of dimensions D simplifies 
considerably in the final expression and only a, D — 1 factor remains. 

We can now combine this information with the full quantum computation of eq. (1121) . we thus 
have 

{V) = e '^^+V^^<e' n(i+VT^)<e^, hi n > 1. (18) 

Remarkably, once written in terms of the de Sitter entropy, the bound on the volume is 
universal, independent of the number of dimensions! As in D = 4 the probability to produce a 
finite volume violating the bound above is super-exponentially small (see [IA\ for details), i.e. zero 
within the effective field theory regime. 



4 Universality of the Bound with Higher-Derivative Cor- 
rections 

The second test of the bound that we provide is with respect to higher-derivative corrections in the 
Einstein-Hilbert plus inflaton action. Higher derivative terms have multiple effects: they modify 
the Einstein equations, the expression for the energy-momentum tensor, the inflaton equations of 
motion and the formula for the entropy as well. Since we are interested here in slow-roll inflation 
we will assume that independently of the modification induced by higher-derivative terms there 
exists a solution where the metric is approximately de Sitter up to slow-roll corrections and the 
inflaton field rolls slowly, i.e. <C H(p. 

This assumption may appear somewhat too restrictive. There are inflationary models (such 
as DBI [19] and ghost [20] inflation) where the scalar field does not roll slowly even though the 
geometry is approximately de Sitter. However, as shown in [13], these models are very far from 
saturating the volume bound as soon they obey the null energy condition and do not exhibit 
superluminal excitations. On the other hand, ghost inflation is capable of violating the bound 
even at the classical level at the expense of violating the null energy condition. Given that the null 
energy condition and the absence of superluminal excitations are crucial for the success of black 
hole thermodynamics [21] we consider this link as yet another indication that the volume bound 
has a thermodynamical origin. We therefore restrict ourselves to models that are perturbatively 
close to slow roll inflation. 
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The most general action for the graviton plus inflaton system can be written as 

^ci^dpuj Rpupaj ^ a-Rpupa J ••■) 0) ~l~ ^kini^Qpuj Rpupaj ^ a-Rpupa y •■•) 0; ^a*^) •■•) ) (-^9) 

where we have isolated the part of the Lagrangian that does not depend on Va4> {^g = j^^R ~ 
V {(p) + . . .) from the one that contains derivatives of the inflaton (£fcm = ~\{d(t)Y + ...). 

We immediately realize that at leading order in slow-roll parameters we can neglect the contri- 
bution of the terms in involving covariant derivatives. This is justified in the leading slow-roll 
approximation after realizing that all the terms of this form are at least of second order in the 
slow-roll parameters. Indeed all such terms must involve at least two derivatives in order to 
contract the indexes. Upon integration by parts, it is therefore possible to have at least two of 
these derivatives acting on different Riemann tensors. The Riemann tensor in de Sitter space 
is proportional to the metric and so covariantly constant, hence its covariant derivatives in an 
inflationary spacetime are proportional to the slow roll parameters. This means that terms in Cq 
involving covariant derivatives start at second order in the slow roll parameters and can therefore 
be neglected in our approximation. 

For similar reasons at leading order in the slow-roll expansion we will only consider up to 
2-derivative terms acting on the inflaton, i.e. 

^kinigpu, Rpupa, 0, Va0, ...) = -^{do,4>dfi(j))U°''^ [q^^, Rf,upa, 4>) ■ (20) 

In the same approximation, since (i90)^ is already subleading in the slow- roll expansion, its coef- 
ficient can be taken at 0-th order in slow-roll, i.e. computed with de Sitter metric, this gives 

n°^(^7^., Rp^pa, ...) = ^7"^n(^7^„ R^,p„, ...) (21) 

where H = Wg^jA. 

We can now proceed to write the Einstein equations. It is possible to show that ignoring higher 
derivatives the gravitation Lagrangian can be written simply in terms of R^'^ pa and R^u^a without 
any explicit dependence on the metric: 

Cg{9pu, Rpupa, 0) = CciR'^'^pa, R^u'a, 0) , (22) 

and analogously for 11: 

l^igpu, Rpupa, 0) = niR^'pa, R^u'a, <P) ■ (23) 

This is proven in appendix |X1 In this way the equations of motion for the metric read 

6C 5R^%^ 5C dR'^fs"^ 1 ~ 1: 



Using the relation 



5C _ 5C 



--g^X = -Ild^<pd,<p. (24) 



SRpupa = (2VpVJg,p + R\pjgpr) , (25) 



^^pvpa '^J^pupa 
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oc=l a=0 



Figure 1: The stretched horizon (in blue) in de Sitter space, with the unit vectors oc and and 
the region T (in red) spanned by the area section T,{t), as defined in the text. 



it is possible to rewrite the equations of motion as (see appendix lAl for details): 

- 2V.V;3— ^ + - Ig^^C = lud^<pd''<P (26) 

where indices in between brackets are symmetrized. 

From here we can see that higher derivative terms change non trivially several relations used for 
the proof of the bound. The Friedmann equation H = — 47rG0^ used to relate the change in entropy 
with the fl parameter gets modified, in particular both the l.h.s. and the r.h.s. receive corrections. 
The modification of the infiaton kinetic term changes the 2-point function of the infiaton in the 
de Sitter phase, which corresponds to a change in the expression of Q. Finally also the expression 
of the entropy changes, since for a generic gravity Lagrangian the Wald formula [22l |23] must be 
employed: 

S = -— [ dJ^^^Q^" (27) 

where k is the surface gravity on the horizon H, dL^y = ^ef^^dA is the area element {e^u is a tensor 
binormal to "H normalized such that e^^^e^" = —2) and 

Q"" = V,e. - 2V,-^e. , (28) 

where is the Killing vector, which on the horizon satisfies the relation V^^;^ = ne^j^. 

We can now proceed to calculate the variation of the entropy in analogy to what has been 
done in the black hole case [21] or in slow-roll for Einstein gravity [H]. We start defining (see 
fig. [1]) a stretched horizon G inside the true de Sitter horizon "H (in FRW coordinates it would 
correspond to the region of points with r = aH~^e~^^ , with a < 1). Within B the Killing vector 
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^'^ = (1, —Hr, 0, 0) has constant norm ^^^'^ = — (1 — a^). The norm is zero on the horizon, since 
^'^ becomes null, and normalized to -1 on the origin r = a = 0, where the FRW observer sits. 
is also characterized by two unit vectors = ^/^/(l — a^)^/^ parallel to the Killing vector and 
= {a, —Hr/a)/{l — a^)^/^, which is orthogonal to B and pointing away from the horizon. In 
the limit where approaches a — )■ 1 and both and become proportional to which 
becomes null. On the stretched horizon, at each moment in time t, we can use these two vectors 



to define the area element (Hl^i, = \{n^Uy — nyU^)dA = ^e^^dA for the 2-sphere which is the 
constant-t section of G. Moving along S spans a three- volume T where k is constant and the 
variation of the entropy can be computed using Stokes theorem: 

An r 4:71 r An r 

- 6— / dE^^Q'^" = f dE^^Q'^" = — / drdAn^V.Q^r (29) 

Js(t) ^ Jar ^ Jr 

where dr is the unit (proper) time interval in B. 

After few algebraic manipulations (see appendix E]) we arrive at the following formula: 

f = ^ / dAU. f 7^^^"/37 - 2V.V,-^) . (30) 
After using the equations of motion f l2B]) we have 



Hdt K H '^^ ^ 3 

The first non trivial result is that in eq. fj3T|) the non-trivial dependence on the Riemann tensor 
of the l.h.s. of the Einstein equation f l24l) is completely gone after using the Wald formula for the 
entropy. We will now show that the residual dependence in the r.h.s. will also disappeared after 
taking into account the modification of the infiaton action. Indeed the kinetic term for the infiaton 
in eq. ( 120|) is not canonical anymore. This implies that the corresponding 2-points function during 
the de Sitter phase now reads 

The corresponding definition (flOl) for Q will thus be 

fi-il^-n2^il (33) 

which finally allows us to write eq. (13T]) again in the form 

f^l2n. (34) 

The differential equation governing quantum fiuctuations of the infiaton will also have the same 
form as in the Einstein gravity case once Q is properly defined according to eq. ( l33l) . 

We thus conclude that, in the approximation we are working, i.e. at leading order in the 
slow-roll parameter, higher derivative corrections in the Lagrangian do not modify the calculation 
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Figure 2: Depending on the starting point, inflaton trajectories infield space (in magenta) drift towards 
different points of the reheating surface (in red), which in general correspond to different vacua. The field 
space is often confined by boundary regions ( in blue ) where, for instance, the energy density of the scalar 
potential become Planckian. 

of the probability distribution of the volume of the reheating surface after inflation. In particular 
the bound on the volume of the universe after slow-roll inflation is universal also with respect to 
higher derivative corrections — the coefficient "1/2" in the exponent of eq. (ITSll does not receive 
corrections from higher- derivative terms! 



5 Universality of the Bound in Multifield Infiation 

After having seen that the bound on the finite volume holds in any number of dimensions and 
after including higher derivative terms in the Eistein-Hilbert action, we now pass to the study of 
the case where we have more than one light field during infiation: multifield infiation. 

It is straightforward to generalize the procedure of ref. [11] that led to eq. ([2]) to the case of 
multifield infiation. We now have 

VV(t, z) - 2VQ, ■ Vf{T, z) + /(r, z) log[/(T, ^)] = , (35) 

where 

V = 9., r = 2vrv/6|, = ^ (36) 

the boldface font is used for vectors in field space, and we assumed for the moment that the 
reheating surface is just a single point at r = r^. We assume there is also a barrier B that bounds 
the moduli space. In this case the generating function must satisfy the boundary conditions 

f{Tr,z)=e-\ (37) 
n .V/(T,z)|,,g = 0. 
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where h is the normal versor to the barrier hypersurface. This is a simple generalization of the 
analogous boundary condition for single field inflation eq. (jlj). The probability distribution for 
the volume is given by 

piV,T)= / dze^'^fir^z) . (38) 

Jo+-ioo 

More generally the reheating region may be a surface in field space TZ, as depicted in fig. [2l In this 
case different points on TZ might correspond to local universes with different physical properties. 
Hubble patches terminating inflation on different points of the reheating region may have different 
Hubble constants, entropy, and even correspond to different vacua (e.g. if 71 is disconnected or 
it represents a moduli space). In this situation it is not clear the meaning of the bound on the 
volume relative to the entropy: the entropy may vary substantially from point to point in TZ, and 
it is not obvious which of the entropies associated to the various vacua we should take. We may 
consider the largest one as an upper bound, but, as we are now going to see, we can have a more 
stringent definition. 

In fact, it is possible to keep track of the type of volume produced at the end of inflation by 
tracking the point where the fields exit inflation. If we assume that the reheating surface TZ is 
made of a set of n disjoint points (r*) in field space, the boundary condition at 7Z for / would 
read 

/(r;,f) = e-^^ (39) 

generalizing the boundary condition of eq. f l37|) by using a different variable Zi for each reheating 
point T^. The conjugate variable to each Zi would be Vi, the volume of type t^. The probabihty 
distribution of creating the volumes V = {Vi, . . . , Vn} of type tv = {t^, . . . , r"} respectively, is 
then the n-dimensional Laplace anti-transform 

p{V,T) = J^dze'^-'fiT^z) , (40) 

where C = {z : z E C A Re(zj) = 0"'"} simply generalizes the contour of eq. ([1]). 

In the case of a continuos reheating surface in field space the vector tv becomes a continuous 
variable G 7Z, £ and V become functions of Tr respectively z{Tr) and V{Tr), and the formula 
for the probability distribution of the volume becomes a functional integral 

p{V{t,),t) = j^T^^Tr) ei^'<'"^<>^<^ f{T,z{T,)) , (41) 

with Tr spanning TZ and normalized so that j^dTr = 1, and f{T,z{Tr)) now satisfying 

VV(t, ^(t,)) - 2 • V/(r, z{Tr)) + /(r, ^(r,,)) log[/(T, z(t.))] = , (42) 

/(r„z(T,)) = e-^(^^), 
n .Vf{T,z{Tr))l^^ = . 

It is possible at this point to show that eq. ( l35l) does not apply only to the case where we 
have a single reheating point. It is also the generating function for the probability distribution 
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of the total volume V = Jj^dTr V{Tr) produced at the end of inflation. Indeed the probability 
distribution for the total volume V is given by 0: 



VV{Tr) p{V{Tr), t) 6 i ¥{4) -V)= dz 6^"" f {t , Z 

\Jn / Jo+-ioo 



(43) 



This is exactly eq. ( I38l) . Therefore we conclude that solving equation (]35l) with boundary condi- 
tions fl371) corresponds to compute the generating function for the probability distribution of the 
total reheated volume, independently of the particular kind of volumes this is made of. 

At this point the equation for the moments of the distribution can be simply generalized 
to the multifield follows 



(nr;)nr;')...nr("^)) = (-l) 



5"/(t,^(t,)) 



6z{t;)6z{t;') . . . 6z{t, 



(44) 



z{Tr) = 



The differential equations for the moments of the distributions are linear and in particular the 
one for the average volume reads 

V'{V{t'^)) - 2Vn ■ V{V{t'^)) + (V(r:)) = (45) 

(nr;))u. = 5(r;-rr), 

n •V(nr;))Ui? = 0- 

The bound: 

We are now ready to give a more stringent definition of the bound that we believe it will still be 
satisfied in general. In the case of multifield inflation, there are multiple kind of vacua. According 
to the meaning of the bound on the volume of the reheating surface in single field inflation, we 
expect that, in the non eternal inflation phase the total volume produced of any kind be bounded 
by the corresponding de-Sitter entropy; in other words we expect a bound to exist for any particular 
kind of volume. We are therefore led to conjecture the following simple generalization of the bound 
on the volume of inflation: 



p(^j dTr V{Tr) > Sup^ [e^^^' ^/^j ; V < +00^ = 

W(r;)p(\/(r;,T)) <Sup^ 



(46) 



r+00 

dTr 

X ieS(-^'-)/25(T/.-T,.) 



^The actual derivation reads 



= r dX [ VV{Tr) I Vz{rr)ei^^<^^^>'^<'^+'''''^i^^<''^<^-^'^ /(^,^(^r)) 

J —OG J 

= I d\ I Vz{Tr)e 

J —00 J 



-iXV 



f{T, z{Tr)) 6{X - tz{Tr)) = -^ / dX e-'^^/(T, ~iX) = / dz e^^/(x, z) 

00 Jo+—ioo 
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Figure 3: Average volume distribution {V{yr)) as a function of yr for the two inflationary examples: the 
waterfall (left) and the tilted waterfall (right). The red-dashed lines refer to the classical evolution. Near 
the classical limit, for large 17, the distribution (in blue) is peaked around the classical exit point. Near 
the phase transition to eternal inflation, at 0,x ^ 1? the distribution (in magenta) broadens (left) and 
drifts towards smaller values of yr (right). 



for any subset X oilZ and with k being a numerical factor of order one. The bound on the volume 
states that the probability to create a volume of kind larger than e^'^^^'^ and total volume finite, 
vanishes up to non-perturbatively small quantum gravity effects. 

In particular, for the average of the volume this implies 

j^dr. (\/(T.))<Sup^[e^(-)/^] , (47) 

and similarly for higher moments. 

Given a certain multifiled inflationary model, one can solve eq. ( H2|) and obtain the probabil- 
ity distribution for any kind of volume after performing the (functional) inverse Laplace trans- 
form (HTj) . In practice, such a task is extremely difficult from a technical point of view. What we 
will do next is presenting two simple examples where we are able to compute exphcitly {V{Tr)) 
and that offer a non-trivial check for the bound in ( H7|) . In fact checking for the average is expected 
to be enough because as it has been shown in the single field case the probability distribution is 
always sharply peaked around the average value. 

5.1 First example: the waterfall 

The first and simplest example we consider is the two-field case r = {x,y) where the slope of 
the potential is uniform everywhere and orthogonal to the reheating surface at = {0,yr), i.e. 
\/r2 = (v^, 0). Physically this corresponds to the infiationary model where the second field acts 
just as a classically irrelevant spectator (Fig. |3]). In this case the two fields decouple because of 
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the shift symmetry in y, and we expect the result to be quite close to the single field case. In 
particular, we expect to recover the same result as in single field infiation once we integrate over 
the position of the reheating point. Though quite similar to the single field case, this example is 
simple enough to allow us to calculate explicitly also the distribution of "different volumes" on 
the reheating surface. This will lead us to quite non-trivial results. 

For example let us calculate the average volume of type yr given the starting point r = {x,y). 
Applying (jHj) we obtain 



6z{yr 



= -^{x,Ay). (4^ 



where we have used the shift symmetry in y to assume that the y-dependence can be only in the 
form Ay = y ~ y^. ip{x, Ay) satisfies 



i^xx + ^yy - 2^Vt^lj^ + V = , (49) 

V^(0,Ay) = -5(Ay) , 
^x{x,Ay\^^^ = , 



where Xb is the y-independent location of the barrier. 
The solution reads 



oo 



^Pix,Ay) = / :^e'>'^yZ±zZ (50) 



— oo 



For f2 > 1 the limit — ?■ oo can easily be done and we get the following expression for the 
average volume of type yr 



{V{yr)) = r 1^6^'=^^+-- = ^^^^^e^^K,[^iQ - l){x- + y-)] , (51) 
J-oc vr + y^ 

where Ki[x] = \ (ite~ 2 (*+t) is the modified Bessel function. 
It is easy to check that the average of the total volume 

/oo /-oo ;-oo ji 2 

dy,(y{y,))= dyr e'kAy+.,_x ^ ^.0-^ ^ ^ ^ (52) 

-00 •/ — 00 J —CO 

coincides with the one- field case, eq. From this we see that the phase transition to eternal 
infiation happens at 

n = l . (53) 

This result is not surprising given the shift symmetry in y. However from eq. flSTj) we also have 
the information on the shape of the volume distribution on the reheating surface. As expected the 
distribution is peaked around yr = y, which is the classical exit point. Let us look at the shape in 
several limits (see also fig. [3]). 



15 



We can expand the Bessel function to obtain 



{V{yr)) e " i+Vi+AyV-^ , for {Vt - l){x^ + /\y') -> I , (54) 

which is a Gaussian in Ay for Ay < a; (when the y-distance of the starting point from the reheating 
point is smaller than the classical trajectory). The tail of the Gaussian turns into an exponential 
for Ay > X while the width of the Gaussian is always smaller than x. 

Instead, very near to the transition to eternal inflation, 1] — )■ 1, the distribution approaches a 
Lorentzian for Ay small enough: 

^^^Vr)) ^ 2^2 ^ (^]-l)(x2 + Al/2)«l (55) 

x'^ + Ay^ 

while the tail becomes again exponential e~^^^~^ ^'^^ for Ay ^ i^t — 1)"^^^. 

The bound on the volume is in this case trivially satisfied. The reheating entropy is the same 
for every reheating point {S{y) = S), and for the volume associated to each reheating region X we 
have 

J dyr {V{yr)) < {V) = e'^^^i+Vi-i/^ < e"^/' . (56) 
5.2 Second example: tilted waterfall 

We now consider a generalization of the former case. We imagine that the reheating surface is 
defined at x = as in the former case. However, now the gradient of the infiaton potential is not 
orthogonal to the reheating surface: VTi = {y/Ilx, V^y)- Consequently the value of the Hubble 
radius at each point of the reheating surface can be different, this fact will allow us to test the 
generalized bound in f H7|) in a non-trivial way. 

In this case eq. (149 p becomes 

'Ipxx +1pyy -2 VUxi^x + Vftylpy + = , (57) 

with the same boundary condition^. The solution for the average volume is similar to the previous 
case: 



, / , - e^^^^^y^yK,[^{n-l){x^ + Ay^)] , (59) 

^yx^ + At/2 



•^Actually since in this case there is a non vanishing slope also in the y direction, the barrier will be naturally 
at angle, normal to the vector vTl, if, for example, corresponds to the region where H =const. In particular it 
will always intersect the reheating surface in one point, as H will be Planckian somewhere on TZ. For this analysis 
we are interested only in the region far from the barrier, where boundary effects can be neglected. So we took the 
barrier to be parallel to the reheating surface as in the previous example and then took the limit — ^ oo. For our 
purposes this is a safe approximation because, as explained in |T3] , finite boundary effects are expected to make the 
actual bound even stronger since they cut-off infiaton trajectories getting arbitrary far from the reheating surface. 
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where we defined 



(60) 



From the solution above, we notice that this time the phase transition happens when becomes 
smaller than one, i.e. before fl reaches one. Indeed if we consider the total average volume: 



(V) 



dyr{V{y)) 



dyr 



dk 
2^' 



(61) 



we can see that it is not analytic at = 1, signaling the onset of the eternal inflation regime. This 
was indeed expected from the symmetries of the problem — it is only the gradient of the potential 
normal to the reheating surface a/^^,. that matters in determining the transition to eternal inflation. 
The slope \/Q.y along the y-direction, parallel to the reheating surface obviously plays no role for 
the phase transition. Still the information about the reheating point is non-trivial in this model, 
as for example each different volume yr is associated to a different Hubble scale or more generally 
to a different kind of vacuum. If we look at the average value of the reheating point defined as 



{Vr) 



IZo dyr {V{yr)) 



y-x 



a. 



(62) 



cl 



we can see that for il^ ^ 1, {yr) coincides with the exit point of the classical trajectory yl 
y — x^JVLy/Vtr^. In this limit, {V{yr)) is a Gaussian sharply peaked around yf^ as in the previous 
example. Instead as VLx — ^ 1, (z/r) — ^ — oo, and {V{iir)) broadens as roughly shown in Fig. [31 This 
result is qualitatively expected: the closer we get to the phase transition, the longer the stochastic 
trajectories become, and the more they drift downhill (because of the \/Vty tilt). This explains 
why the singularity in (V{yr)) only appears at f2 = 1. Indeed even when VL^. < 1, {Viyr)) can 
be finite for any y^., the infinity of (y) being due to trajectories exiting at y^ — ?■ — oo. When the 
whole f2 = 1, then also (V^(?/r)) has to start diverging because trajectories can start to go up-hill 
also in the y direction and produce infinite volume at finite yr- Graphically the phase diagram 
looks like as follows: 



(y)=oo 

{V(y))<cc 



{V)=oo 




(y)<oo 
(y(v)><oo 



> 
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Notice that for < ^ but > 1 (the yellow region above) we have eternal inflation only because 
the total volume diverges, while {V{yr)) stays finite for every finite yr- In realistic models we 
would expect also the range of y to be finite (as the shift symmetry in y is broken by Qy). In this 
case the divergence in (V) disappears being due to the non integrability of {V{yr))- The yellow 
region turns into a non-eternal slow-roll inflation region and the phase transition is expected to 
happen now at Q = 1. 

We are now ready to formulate and test the bound in this multi field case. Eq. ( 1471) must be 
checked for any interval X on the reheating surface, however since the entropy in this example is 
monotonically decreasing with yr it will be enough to consider just the intervals I = [y^, oo) for 
every y^. Notice that the total volume is bounded by 



{V) = e(v^--v^)- = e ^+V^ < (63) 

where Sc is the entropy as calculated at the value of the Hubble constant corresponding to the 
classical exit point y^^ of the reheating surface and we used the relation x = Gy/UxNc- It follows 
that the bound is trivially satisfied on any interval X = [?/°, oo) with y^ < yc since 

dyr {V{yr)) < {V) < < ^s(v^r)/2 _ (54) 



We just need to show now that the bound is not violated when y^ > yf. 
For {VL - l)(a;2 + Ay^) > 1, (F(?/r)) behaves like 



which is a falling exponential for > yf (since (^(y-r)) is peaked around (y.,.) < yf). Therefore 
the integral in the l.h.s. of eq. fH7|) is dominated by the value of the integrand at y^- We thus 
need to check whether 

{V{yl)) < , (66) 

up to pre-exponential factors. At leading order in the slow roll approximation we also have that 
^hjr) = Sc — 2^J^Ly{y^. — yf) and looking only at the exponents of fl66l) we have 

^/n'^x + v^Ay - v^(?T^1)(^^Ta7) <^- v^(y° - yf) . (67) 

Using the fact that x = Gy/Tl^Nc and defining Sy = y^ — yf, after a bit of algebra we get 



This is the sum of two terms. The first is positive and the second is also positive since Sc > AS 
Sc ~ Sstart = 12fiA^c- 
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Of course, the bound continues to hold even when {Q — l){x'^ + Ay'^) 3> 1 is not satisfied — right 
before entering the phase of eternal infiation — because the distribution of volumes move further 
towards large negative values of i/r where the bound is trivially satisfied because the volume 
produced in X is smaller. 

Notice also that the l.h.s. of ( 168|1 is minimized at 6y = Aye (corresponding to = y), where 
we get 



VL^\/VL - 1 + n — ^ 1 



> , (69) 



which can be saturated only when Sc = AS and Q = 1, which implies Qy = 0, i.e. the untilted 
potential of the previous section, corresponding in practice to single-field infiation. 

This shows that tilting the reheating surface with respect to the gradient in field space actually 
makes the bound stronger. It also suggests that only single field slow-roll infiation manages to 
actually saturate the bound, any other multifield models which move further away from the 
symmetric setup of the single field case, seem to make the bound stronger, by producing less 
volume. 
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A Some explicit computations 

In this appendix we report some of the calculations omitted in section |H 

First we show that any given Lagrangian C, generic function of the metric g^^, the Riemann 
tensor R^upa and scalar quantities 0, with no explicit covariant derivatives, can be rewritten just 
in terms of the Riemann tensors with two covariant and two contravariant indices [R^'^ and 
R^u''u) without any explicit dependence on the metric g^u, i-e. 

C{gp.. Rp^pa, 0) = £(i?^V' <P) ■ (70) 

Any term in C can be viewed as a network of nodes and lines, each node corresponding to the 
insertion of a Riemann tensor, each line to the contraction of two indices with a metric. Each 
node has thus four lines attached to it. The statement f lTU]) corresponds to showing that the lines 
of each network can be oriented such that each node has two incoming and two outgoing lines (see 
fig-i). 
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Figure 4: Each node corresponds to the insertion of a Riemann tensor, each link to the contraction of 
a pair of indices. Closed orientable loops can be identified such that all nodes will end up having an 
equal number of incoming and outgoing lines. Hence any scalar contraction of Riemann tensors can be 
written just in terms of Riemann tensors with two upper and two lower indices contracted without using 
the metric. 



Consider now a node of the network and let start following a line departing from it. Since 
every node has an even number of lines attached, when the line arrives to a node there will be 
another line which has not been used yet, which can be followed, until eventually the path will 
close arriving back to the initial node. At this point we have a closed path which we can orient. 
Each node in this path will have an equal number of incoming and outgoing lines and an even 
number of unused lines. We can thus repeat the construction by starting with another unused 
line, forming another closed loop and orienting it. When we will have used all the available links 
we will have oriented all lines such that an equal number of incoming and outgoing line will pass 
through each node. This proves eq. ( TfOl) . 

We prove now that the equations of motion can be rewritten as in eq. (l26ll . First we use 
that 

£(i?^V, ^V.,...,0) 
6C 



6C 

where Sn/SnRf_iupa means that we only differentiate with respect to the R,_iupa appearing in the n-th 



6iC 



^iRafipcr 

hi 

^iRoLvRa 



9ap9pv 
9ap,gi3p 
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argument. We thus have that 



dK"Ppa OK'^pl'a OlKaPpa d2KauPa 



^l-Rpupa ^2-R'pupcr 

6C 



^Rfii/pa 



{2VpVJg,p - R%pJ9ap) (71) 



where in the last step we used the identity (125|1 and the fact that 



^lRp,upa ^2Rpupa ^Rfiupa 

Hence, we finally have 

6C 5R-^p. 5C 5R%''. _ SC 6C 

5R-^. 5gp. 5R-pP^ 5g,, " " 5Rp^,,)„ ' ^ ^ 

which can be used to derive eq. fj26|) from eq. fl24|) . 

The last missing step to explicitly show the universality of the bound with respect to higher 
derivative terms, is the relation between the variation of the de Sitter entropy and the metric 
equations of motion, in particular how to get eq. ( 130|) from the variation of eq. (l27|l . We start 
with the Wald formula for the entropy of ^(t), which reads 



Sit) = / dS^, — V,^. - 2V,— , (74) 

pupa '^^^pvpa / 

where k is the surface gravity on the stretched horizon. Using Stokes theorem it follows that 

dS__4n r d^^^^ ^ f-^^V e - 2V ^'^ ^ 

dt K Jx dt ^ " \5Rpyp„ ^ " ^ ^Rpvpa 



An 



dA np_ 



'^^pupa '^J^ppua / '^^pupa 'J-^pupa 



- I dAh.i^ ( -^R\pr - 2V.V 



T \'-'J^pupT '-'J^pupa 

I dAU, (j^R'^pr - 2V.V,-^ ) , (75) 

where in the first step we defined = = {a, —Hr /a, 0,0), which is normalized like the 

Killing vector [n'^np = —C.'^C.p), in the second step we used the cyclic properties of the Riemann 
tensor, in the third the fact that 

V(^^.) = , V^V^^p = R%.p^. , (76) 

and in the last we performed the horizon limit a — ?■ 1. Eq. fITS]) matches eq. fl3UI) and this 
terminates our proof. 
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